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Thermodynamics and kinetics of lithium intercalation into C-SiC>2 nanocomposites are investigated. Depen- 
dencies of both differential capacity and intercalation kinetics on the nanocomposite size are established. The 
processes are analyzed in terms of the impedance model. The obtained results are explained based on the 
quantum effect of interference blockade of electron tunneling into a nonmetallic nanoparticle. Propositions for 
the new electrochemical energy storage technology are presented. 
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1. Introduction 

A growing interest to nanostructures has been observed lately. These structures possess a num- 
ber of unique physical properties that turned out to be quite promising from the viewpoint of their 
practical applications in electronics. For instance, nanotechnology may be one of the effective ways 
of solving an urgent technical problem of producing a cathode material with high specific energy. 
Really, application of the nanodispersed FeS2 in an energy storage device with the lithium anode 
increases the specific capacity by about 20% in comparison with the coarse-grained homologue [T] , 
and the nanosized a-Fe 2 03 possesses high recirculated capacity 200 mA-hour/g and good cycling 
in the range of 1.5-4.0 V regarding Li + /Li in comparison with the macrostructured a-Fe2C>3 , 
a-Fe 3 4 , and 7-Fe 2 3 0. 

On the other hand, the urgency of searching for new cathode materials is closely connected 
with lack of raw materials for traditional cathode-active ones |3J . We think that production of new 
cathode materials on the basis of cheap and ecologically clean substances due to the dimensional 
(nano) effects is a rational way of making the energy storage devices. Beneficial application thereof 
obviously depends on the degree of understanding the physical processes in nanostructures. The 
paper is devoted to consideration of the above problem. 

2. Samples and experimental setup 

Let us consider a silicon dioxide Si02 which is a widespread ecologically pure material. The 
nanosize silicon dioxide intercalated by Li + belongs to the groups of materials usable in developing 
the high capacitive energy storage devices. 

It is well known [3] that silicon dioxide occurs in three structural forms, i.e., quartz, tridymite 
and cristobalite. None of them in macroscale is usable as a cathode-active material [5]. This is 
connected with the extremely low level of density of ionised defects at a room temperature and 
low concentration of the shallow level of trapping, although "guest" positions for the introduced 
lithium are well defined by nanodispersed tridymite structural channels (figure [TJ . The situation 
has drastically changed with the transition to the nanosize SiC>2 with lithium introduced into its 
structural channels. In this case, the density of states at Fermi level can become considerable due 
to surface states (the "electrochemical grafting" concept |5J|3]) and due to sharp decrease of the 
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diffusion resistance. Our quantum-mechanical calculations [7 have shown that the lithium pene- 
tration into the nanodispersed tridymite structural channels (figure [T]) is an exothermal process. 




Figure 1. Structural channel in the tridymite cluster along [0001] (parallel to C-axis) (a) and 
along [1110] (normal to C-axis) (b). 




The nanodispersed silicon dioxide used throughout experiments was obtained from the sili- 
con tetrachloride by pyrogenic method in hydrogen-air jet in the reaction system with minimum 
turbulence at 1000h-1500°C. 

The size of the sedimentation fractions was an- 
~ . ~ PT"iij"i alyzed both by small-angle X-ray dispersion [5] 

'-- w and by electron microscopy. Electrodes with sur- 

face area 0.5 cm 2 on a nickel substrate were formed 
for electrochemical study. Composition of an elec- 
trode is defined by ratio active material (SiC^): 
conductive agent: cementing agent as 85% : 10% 
: 5%. The silicon dioxide mass does not exceed 
1 mg. Structure of the electrodes is schematically 
shown in figure [2] 

Figure 2. Fragment of the electrode structure Thermodynamics and kinetics of the lithium 

(TSG - thermosplitting graphite). intercalation were investigated in the three- 

electrode electrochemical cell with the monomolar 
LiBF4 solution in 7-butyrolactone. The impedance analysis was performed within the frequency 
range 10~ 2 -^10 5 Hz using AUTOLAB device manufactured by ECO CHEMIE (Holland) equipped 
with FRA-2 and GPES computer programs. Galvanostatic "charge-discharge" cycles were realized 
by standard electronic circuit. 

The electromotive force method [H] was used for thermodynamic analysis of the process. Figure[3] 
presents the dependence of differential capacity dx/dU rc i on a; (a; is a number lithium per formular 
unit of a host material; U le \ is a relaxed value of voltage of the open circuit relative to lithium 
electrode depending x content). It is seen that the current-generating reaction only for the nearly 
9-nm nanoparticles undergoes the I-st type phase transition (jump of dx/d[/ ro i) with the formation 
of a two-phase state. At x defined, X-ray investigation fixes in LiSi02 the NaCl-type phase in 
which oxygen anions form a face-centered cube and the ions Li + and Si 4+ occupy the octahedral 
hollows. The lattice parameter of a such phase is 4.0287±0.0009 A. Figure [3] shows a researched 
current-generating reaction in the nanodispersed silicon dioxide [7] with the average size of the 
particles of the order of 16, 11, 9 and 5 nm intercalated with lithium content of 0.8 < x < 2. It is 
seen that the specific discharge capacity, specific energy, at x — x max for which discharge voltage 
is 1.5 V, essentially depend on the nanoparticle size (table [l}. The maximum flow of the reaction 
takes place in case of the 9 nm size. What exactly has caused the situation? 

In order to analyse the Faraday's current-generating reaction in various nanosized silicon diox- 
ide, let us consider Nyquist diagrams presented in figure [4] They show a strong increase of the 
resistance R c of the charge transport (figure [4] curves 2, 4) for all sizes of nanoparticles (with 
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Figure 3. Gibbs function dependence on the guest lithium concentration in the silicon dioxide 
nanosystem with average size of particles 16 nm (1), 11 nm (2), 9 nm (3), 7 nm (4) and 5 nm (5). 

Table 1 . Specific discharge capacity. 
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Figure 4. Nyquist diagram of the Li + -intercalated current-generation for the silicon dioxide 
nanosystem with average size of particles 16 nm (1), 11 nm (2), 9 nm (3), and 5 nm (4). 
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the exception of the 9-nm size in which the kinetic mechanism takes place; curve 3) up to the 
practically full blockade of the process (curve 1). In other words, the interphase charge transport 
is a well-marked effect. 

Herein below we present particular considerations that will enable us to comprehend the cause 
and circumstances at which this situation becomes possible. 



3. Model 

Faraday's current-generating reaction is connected with penetration of an electron into the 
nanoporous hollow filled with lithium. Today there are some investigations of the behavior of elec- 
trons in nanostructures. The researchers frequently use well-known quantum mechanical models. 
In particular, the effect of a resonant laser impulse irradiation on an electron within the system of 
two semiconductor dots was researched in [10]. Two interlinked rectangular wells that model the 
quantum dots were used. The paper shows that quantum dynamics of electron under the action of 
electromagnetic irradiation of nanostructure essentially depends not only on the behavior of the 
wave but also on geometric and energy characteristics of the dots. Resonant transition in three- 
barrier nanostructure was studied in |11| . Therein a one-dimensional model was used with three 
rectangular wells of various depth, i.e., like in paper [10] . the potential on the whole is a nonanalytic 
function. However, the solutions of the problem (wave functions, energy levels) are known for each 
well. The effect of rearrangement of levels in the system may be taken into consideration from the 
procedure of matching the wave functions together at the border discontinuity of the potential. 

As it follows from the above cited papers, as well as from several others, the present state of 
theoretical research of nanoobjects justifies the application of even the simplest models for the 
purpose of comprehending the physical effects therein. 

Below we make a research within the simplest model of energy storage in nanosystems, namely, 
the electron tunneling into a nonmetallic nanoparticle in a) one-dimensional model "border-core- 
border" and b) model of porous nanostructure (or two-dimensional nanoparticle, nanotube) 



3.1. One-dimensional model "border-core-border" 

In figure [5] within the framework of one-dimensional problem, the nanoparticle is presented by 
potential C/ 2 of a hollow with two potential barriers U\ and U 3 which set off the hollow. Different 
thickness of such potentials should be regarded as effective potentials of particles of different sizes. 
Then electron with energy U2 < E < min{J7i , L^}, moving from left to right, will be tunnels 
through a barrier a into the particle. The barrier c shows itself mediately, namely through the 
electron de Broglie wave reflected from it. 




a b c 



Figure 5. Onc-dimensional potential simulated a nanoparticle. 



It is well known that at certain conditions, the tunneling probability of an electron through the 
system of two identical rectangular potential barriers reaches unity [12| . This is just the essence 
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of resonant tunneling. This situation is observed for a subbarrier electron with the energy value 
coinciding with its virtual values in a well. This is caused by cancelation at superposition of the 
tunneled electron de Broglie wave into a well and the elastically reflected wave from the opposite 
barrier. Contrary to the traditional resonant tunneling problem through the both barriers, our 
interest was focused on tunneling into the hollow, i.e., tunneling through one of them, which turns 
out to be the essence of current-generating reaction. 

Due to non-analyticity of the potential, it is presented as the sum of potentials: 



0, x e (-oo,0], (0) 

U u xe[0,Xi], (1) 

U (x) = { U 2 , x G [xi ,x 2 ] , (2) 

Us, x G [x 2 ,x 3 ] , (3) 

0, xe[x 3 ,oo), (4) 



(3.1) 



i.e. all area is broken into intervals, in which the potentials are analytical functions (a = xi 
b = x 2 - x\ ; c = x 3 - x 2 ). 

The solutions of the time-independent Schrodinger equation are of the form: 



ip s (x) = C s exp(k s x) + b s exp(—k s x) 



(3.2) 



with k s = \\/2m (E — U s )/H 2 ; index s 
figure [5] 



- 0, 1, 



fen 



, 4 coincides with the notations of intervals in 

const • i ■ \[W- 



k,4 = const 
U 2 , where const : 



'E, k x 



Ui. 



Thus, taking into account potential (3.1l 
k 3 — const • i • \/E — U3 , k 2 = const • i ■ \/E - 

In our case, tunneling through a barrier a is described by the tunneling probability D = 
\G 2 jCo\ — C 2 C 2 /CqCq . Thus, it is necessary to find the weight multipliers C 2 , Cq of the corre- 
sponding wave functions. Taking into account the properties of the wave function (the absence of 
both the jump of the function and its derivative in any point) and matching such functions and 
their derivatives in the jump points of the potential one, yields a set of eight linear equations in 
Co, bo, Ci, 61, C 2l b 2 , C3, 63, C4 (in ip4 (x) 64 = due to the absence of the reflected wave in 
this interval). Analogue of such a system is shown below. 

A standard procedure of numerical solution of the equations was used. In figure [6] (a) the 
logarithm of the tunneling probability on the nanoparticle size is presented at parameters a = 5 nm, 
c = 4 nm, Ui = 4.0 eV, U 2 = 1.0 eV, U 3 = 5 eV, U 2 = 1.0 eV, E = 3.0 eV, and in figure [6] (b) - 
at the same parameters with the only replacement of E = 3.0 eV by E = 3.9 eV. The replacement 
also describes the sub-barrier electron but with energy closer to the lower barrier top. 
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Figure 6. Dependence of the tunneling probability on the nanoparticle size at E = 3.0 eV (a) 
and E = 3.9 eV (b). 
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From figure [6] it follows: 

• The property of the tunneling probability D as a function of the well thickness [In I? = / (&)] 
is its oscillation behavior. In the points of minimum, D is practically equal to zero. Such 
behavior is a result of superposition of the wave tunneled into b and the elastic wave reflected 
from the barrier c. 

The value of the tunneling probability essentially depends on the ratio between the energy 
of electron and the height of barriers. Comparison of figure [6] (a) and figure [6] (b) shows that 
change of not only the value but also of periodicity of the oscillation In D — f (b) . 

• The increased thickness of an opposite barrier (c) decreases the tunneling probability without 
changing its periodicity. 

Thus, depending on the ratio between geometrical and energy characteristics of the model, 
the electron may or may not penetrate into the hollow. The actual situation may be somewhat 
distorted if one takes into consideration the Gaussian-like distribution of nanoparticles by size as 
well as energy distribution of the bombarding electrons. However this does not change the main 
qualitative conclusion of the analysis, i.e., the oscillation dependence \nD = f(b). 

3.2. Porous nanostructures (or two-dimensional nanoparticle, nanotube) 

Consider electron tunneling in porous nanostructures (two-dimensional nanoparticle, microp- 
ores, nanotube in two models: a) with the potential 

f 0, x 2 +y 2 > R 2 , 
U(x,y) = l U u r^^Ty^^R, (3-3) 
[ U 2 , x 2 +y 2 < r 2 . 

Here, space is divided into areas in which the potentials of the electron are analytic functions. 
The model potential is presented in figure [7] Therein, a pore U2 is restricted by a rectangular 
barrier of U\ in the form of a ring with an outer radius R and inner radius r (U2 <U\). 




Figure 7. Two-dimensional potential simulating nanostructure pores. 
In the model b), the potential is as follows: 

0, x 2 + y 2 /10> R 2 , 
U(x,y)= { U u y/x 2 + y 2 /10 ^ R, 

U 2 , x 2 + y 2 /10 < r 2 . 



(3.4) 



Such a model resembles the model of cylindrical tube [model a)] with the sole difference being 
that in this case its cross section is an ellipse. Comparison of (3.3 1 and (3.4 1 shows that the 
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short semi-axes of ellipses, that restricts the hollow along the OX axis, coincides with the radii in 
model a) . Large semi-axes of a size y/WR and y/Wr are oriented along the OY axis. 

We analyze the manifestation of resonant tunneling in both models. Herein, solutions of a 
time-independent Schrodinger equation are of the form: 

i>s(x, y) = C s (y) exp(k s x) + b s (y) exp(-k s x), (3.5) 

Here, as well as above, k s = u/2m (E — U s )/h 2 , while index s = 0, 1, . . . , 4 coincides with the 
notations of the ranges in figure |7| Thus, fen = &4 = const • i • \[E. k x = k 3 = const • i • yJE — U x , 
k 2 = const • i ■ \J E — U 2 , where const = \J 2m / ti 2 . The presence of x and y in the argument of wave 
function shows the two-dimensionality of the problem. 

Let the electron with energy U 2 < E < U x be moving from left to right along the Ox axis. 
Herein, the passage through the barrier is possible only due to tunneling. In the present case, the 
barrier can be considered as a set of one-dimensional barriers for fixed values of y. Geometric sizes of 
the barrier are varied depending on the specific value of y. Hence, for example, in model a) at y = 0, 
the barrier is the thinnest, (R — r), while the size of the hollow is maximum, (2r). In the second 
limiting case, at y — r, two opposite barriers interflow and their overall thickness becomes equal 
to 2yR 2 — r 2 , while the hollow thickness becomes zero. In general, a traveling electron parallel to 
OX, for a fixed y will be reflected from barriers at an angle different from tt. The larger is y the 
smaller is the angle. Below, we restrict the analysis of tunneling in the y area to no more than 
the neighborhood of 0.01% of the pore size along OY relatively to y = 0. (As will be seen, the 
y = neighborhood plays a key role in tunneling in these models) . In this case the reflection angle 
practically equals 7r. 

The tunneling probability of a homogeneous flow of electron through the left-hand barrier 
into a hollow of nanoparticle is D (y) = |C 2 (y)/C (y)\ 2 — C 2 (y) C 2 {y)/Cg (y) C (y), while the 
total probability is determined through the integral by y from up to r of this characteristic. 
Here, Co (y) is the amplitude of an incident wave, while C 2 (y) is the amplitude of the wave that 
has transmitted this barrier. These amplitudes may be derived from the system of eight linear 
equations in respect to C (y), b (y), C x (y), b(y), C 2 (y), b 2 (y), C 3 (y), b 3 (y), C 4 (y) following 
from the continuity condition of wave functions and their derivatives in any point, including the 
points of potential discontinuity. At matching together the wave functions ipo{x,y) and %j)x{x,y) in 
the point xi(y), which is common for the regions and 1 (it is the root with minus sign of the 
equation \/ x\ + y 2 = R) , we shall get the following: 

C (y) exp(fc • xi) + b (y) exp(-fc ■ x x ) = C x (y) exp(fci • x x ) + b x (y) exp(-fci • x x ). (3.6) 

Similarly, matching together the wave functions %j) X {x,y) and i/}2(x,y)in the point x 2 (y) (it is 
the root with minus sign of the equation yj x\ + y 2 — R) yields the equation: 

Ci (y) exp(fci • x 2 ) + b x (y) exp(-k x ■ x 2 ) = C 2 (y) exp(fc 2 ■ x 2 ) + b 2 (y) exp(-fc 2 • x 2 ). (3.7) 

Similarly, for the points x 2 (y) and x 3 (y), and X3 (y) and X4 (y), we shall get corresponding 
equations: 

C 2 (y) cxp(fc 2 • i 3 ) + b 2 (y) exp(-fc 2 • x 3 ) = C 3 (y) exp(k 3 ■ x 3 ) + b 3 (y) exp(-fc 3 • x 3 ), (3.8) 

C 3 (y) exp(fc 3 • x^ + b 3 (y) exp(-fc 3 • x 4 ) = C 4 (y) exp(fc 4 • x A ) (3.9) 

where x 4 = —x x , while x 3 = —x 2 . 

Continuity of the first derivatives in those same points yields still four equations: 

C (y) k exp(k -x x )-b (y) k exp(~k -x x ) = C x (y) k x exp(k x -x x )-b x (y) k x exp(-k x -x x ), (3.10) 

Ci (y) k x exp(ki-x 2 )-bi (y) k x exp(-k x -x 2 ) = C 2 (y) k 2 exp(k 2 -x 2 )-b 2 (y) k 2 exp(-k 2 -x 2 ), (3.11) 
C 2 (y) k 2 exp(k 2 -x 3 )-b 2 (y) k 2 exp(~k 2 -x 3 ) = C 3 (y) k 3 exp(k 3 -x 3 )-b 3 (y) k 3 exp(-k 3 -x 3 ), (3.12) 
C 3 {y) k 3 exp(fc 3 • x 4 ) - b 3 (y) k 3 exp(-fc 3 • ar 4 ) = C 4 (y) fc 4 exp(fc 4 • x A )- (3.13) 
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Within the system of eight homogeneous linear equations ( 3.6 H 3.13 1 there are nine unknown 



quantities. However, having divided all the equations into Co(y) we get a system of non-uniform 
equations with eight unknown quantities, including ^[^y, which is connected with the tunneling 
probability into the hollow. 
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Figure 8. Dependence of the tunneling probability In [D(y) / Do] on y for pores with ring (a) and 
elliptic (b), (c) profiles. 



Calculations of the tunneling probability were carried out using the program Maple-13 with 
the parameters, U\ — 40 eV, U2 = 1 eV, E = 4 eV, R — 10 nm, r = 8 nm, and semiaxis of 



in the form f(y) 



D(y) 



where 



s in the same range of y, and in 



the ellipse b 2 — 10. The obtained results are shown in figure 
ln_D = 51.320595. Figures [8] (a), (b), show f(y) in both model: 
figure [8] (c), the dependence f(y) in the model (b) in a wider range of y. In such y range in the 
model a), f(y) has a similar form but with less marked nonmonotonicities. Due to a large range 
change of f(y), its real form is somewhat garbled by the applied Maple-13 spline. Nevertheless, 
even this form allows us to observe the manifestations of resonant tunneling in the models. 



4. Discussion of results. Conclusions 

The obtained dependencies In = f(y) (see figure |sj) are generally of decreasing character 
with a set of repeating plateaus. 

The decrease of tunneling probability with an y increase is apparent. For any fixed y it corre- 
sponds to the one-dimensional eigenvalue problem with two identical square barriers. The increase 



23705-8 



Quantum mechanic tunneling 



of parameter y in our case causes the widening of the barriers and at the same time the narrowing 
of the hollow. It is known that in case of tunneling through a single barrier, the value of tunneling, 
similarly to our problem, decreases but this dependence does not have plateaus, which we relate 
to the resonant tunneling. 

Let us consider the virtual states of an electron in a hollow, i.e., between two barriers. In order 
to qualitatively comprehend the obtained behavior of the tunneling probability, it is sufficient to 
employ the results of a quantum-mechanical problem of a particle in infinitely deep square well. It 
is well known that here the energy spectrum of a particle is a set of discrete levels E n = |^n 2 
(in is the mass of a particle, n is the principal quantum number, n = 1, 2, 3, . . . ). In our case, 
the dependence of the well width d on y (and consequently E n ), causes the observed behavior of 
the tunneling probability. Compare the behavior of the ground states E\ (y) for two types of the 
pore cross section, i.e. the circular and more general elliptic one extended along OY. In the first 
case, E la (y) = 2m J r2 h _ y 2 ) , and in the second case E lb (y) = 2 m[r^-(v/b) 2 \ ' ( b ^ 1 is tne semi - axis 
of an ellipse). It can be seen that upon the axis, i.e., at y = 0, these levels coincide. However, the 
change y — > y + Sy causes a more rapid rise of position of the level Ei a (y) in comparison with the 
growth of En,(y). Within the limit b — > oo, En,(y) does not depend on y, i.e., this will be a level it 
"degenerates" by y. In a general case, a step-by-step transition y — > y + Sy is accompanied by the 
appearance of closely positioned levels. 

As noted above, at y = 0, the values f(y) coincide for both models. In the present case, f(y) 
is nonzero. It means that the energy of the incident electron E differs from the virtual levels. 
Suppose that E is between the virtual levels E n and E n+ \. Increasing y (which is accompanied by 
narrowing the cavity) the virtual levels rise on the energy scale. In model a) the n-th level E(y) 
with y ~ 0.00029 coincides with energy E of the incident electron. As a result, there is a sharp 
drop f(y), i.e. the resonant tunneling effect. A similar behavior of f(y) is observed in the case b), 
but for E(y) at y ~ 0.00093. The difference between these values y is connected with the fact that 
the hollow thickness in the model a) changes more rapidly with the change of y than in the model 
b). A further increase of y, for example, in the model a) is accompanied by resonant tunneling due 
to the levels E n at E(y) for y up to its value ~ 0.0006. At y ~0.0006, the next virtual level, 
coincides with energy E of the incident electron. Further, the non-monotonous series of /(y) follows 
up to y ~ 0.0008 at which the next virtual level, E n _2, coincides with E and so on. 

The main conclusion of the paper consists in stating that the degree of the electron tunneling 
into the hollow is determined both by the geometry of the pore and by the energy of the bombarding 
electron. 

The obtained results may be represented alternatively. Let us consider a homogeneous flow 
parallel to the axis y = of monochromatic electrons. In the range of y where virtual states 
coincide or are close to the energy of bombarding electrons, their contribution in tunneling into 
hollow is either altogether absent or just minimal. A more real case is the flow of electrons with 
some distribution by energy. Assume that this is Gaussian-type distribution with a maximum at 
E = £ max . Even if E max coincides with virtual level, there are the energies from the "tails" of 
the Gaussian that do not coincide with this level. Exactly these electrons ensure the tunneling. 
However, since the number of such electrons according to the Gaussian distribution is low, their 
contribution into tunneling will be negligible. In other words, a more real picture being taken into 
account does not change the qualitative conclusions obtained earlier. 

Herein above we have focused on the maximum tunneling taking place at y = in case of a 
single barrier. The presence of the second identical potential barrier and hence a possible resonant 
tunneling, somewhat changes this statement. In the case of geometric parameters of the pore, such 
that at y = 0, its virtual levels coincide with the energy of the electron, this electron does not 
provide any contribution into the tunneling. The tunneling is realized by electrons that are more 
distant from y = and their energies do not coincide with the virtual levels. Thus, due to a larger 
thickness of barriers in comparison with the y = case, the degree of tunneling will be much 
lower. However, suffice it to change the energy of the electrons so that it does not coincide with 
the values of virtual states in the region y = 0, then the degree of tunneling will abruptly grow. On 
the other hand, the degree of tunneling may be essentially different for the electrons with the same 
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energy depending on the geometry of the pores. In other words, the degree of tunneling is not much 
effected by the energy of electrons or by the geometry of pores, but rather by their interrelation. 
Thus, the above mentioned result of the paper [7 concerning the maximum tunneling in 9-nm 
pores in comparison with the 5, 11, the 14-nm pores may be caused by the energy bombarding 
electrons which least of all coincides with their virtual states. 

The obtained results make it possible to propose the impedance model of the investigated 
mechanisms presented in figure [9] (a) . Here i?E is the resistance of electrolyte, R A is a relaxation 
resistance of the low-resistance phase A and is its chemical capacity (chemical capacity is the 
ratio of carrier concentration n and Fermi level Ep n , C M = e 2 |14|). Due to quantization of 
an energy spectrum in nanoparticles, the chemical capacity is transformed into the quantum ca- 
pacity [12], and it noticeably effects the transport processes [15] . This effect should be expected 
only in nanoscale objects because in this case only the quantum capacity can disconnect the circuit 
R3-C3 (figure [7] (b)), i.e., the DC-mode. 
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Figure 9. Equivalent electric circuits of impedance model for the Li + -intercalation in the nan- 
odispersed silicon dioxide: a - common, b - at the interference blockade of electron tunneling, 
c - in presence of electron tunneling. 



The branched R A t s || C AS circuit produces a charge transport through the barrier a and 
Rf || Cf - through a barrier c, which divides silicon dioxide and electrolyte. The Rendls-Ershler 
circuit [16] produces the lithium cations transport from electrolyte into the structural channels 
Si0 2 . 

According to the results of theoretical calculations it is obvious that R A t s (as well as C AS ) 
depends on the particle size. Therefore, at maximum and minimum value R A t s , the equivalent 
circuit 9 (a) transforms into the circuits 9 (b) and 9 (c) (see figure |9|. It is easy to see that the circuit 
9 (b) (figure |9| is capacitive. It well describes the impedance behavior at the electron tunneling 
blockade into the nanoparticle Si02 (figure [4] curve 1). The circuit 9c corresponds to deblocking 
the electron tunneling and to the formation of two-phase state. Results of the proposed model well 
agree with a package of experimental data: Kramers-Kronig test does not exceed 3 • 10 -5 , and the 
frequency dependencies of ordinary difference are of completely random character. According to 



the computer parametric identification (figure 10 1, the difference between experimental and model 
curves does not exceed 5%. 

Figure [lO] shows an essential dependence of the parameter R ct on the particle size. Its sharp 



minimum at d = 9 nm coincides with the anomaly 



d(7 ro 



at this value. Here is another interesting 



fact. Intercalated lithium into nanoparticles of various sizes (excluding very small ones) increases 
the charge transport. 

It should be noted that the absence of diffusion control in this case and, as a consequence, 
the minimum content of the diffusion resistance in the total resistance of a system, is extremely 
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Figure 10. Dependence of the parameter R c t on the particle size. 

important to practical application of the 9-nm particles in the current-generating process 

SiC>2 + xLi + + xeT = Li^SiC^ 

with a horizontal discharge curve at 0.9 < x < 2. By contrast, another situation for the 5-nm 
particles (obviously, for the 16-nm particles as well) is observed - the positive charged lithium 
is screened by electrons, injected from contact to a conductive agent. Such an effect (obviously, 
possible only for nanoscale objects) may be defined as sub-Faraday or pseudo-Faraday discharge 
with the expected hypercapacity. Actually, at discharge of a galvanic cell with electrodes on the 
base of both the 5-nm and 16-nm silicon dioxide, the horizontal plots in voltage-time characteristics 
are absent. It testifies to the process without the 1-st type of phase transition. Results of the 
potentiometric measurements do not allow us to claim the prevalence of the diffusion control over 
the kinetic one. At least in the range of 2.94-1.85 V, the dynamic current- voltage characteristic 
is close to ideal polarizability. The calculated specific capacity of "a condenser regime" in a given 
electric voltage range for nanocomposites on the basis of the 5 nm silicon dioxide was ~T075 F/g, 
and for the 16 nm one it was 2034 F/g. Taking into account the active surface area, differential 
capacitance reaches the record value ^256 and 2400 /iF/cm 2 , accordingly. The latter value exceeds 
the known maximum value of the differential pseudo-capacitance by about five times [17| . As in the 
ideal polarizability case, the multifold charge-discharge cycles are possible here. This emphasizes 
the importance of the above considered nanomaterials to the technology of electrochemical storage 
with the high cycling. A similar conclusion is less obvious for the 11-nm nanosize titanium dioxide. 
Obviously, faraday and pseudo- Faraday processes co-exist here. 

5. Conclusions 

1. The nanosize silicon dioxide with the Li + -intercalated current-generating reactions in it is 
suitable for technology of electrical storage with the essentially higher specific capacities and 
energy than in the traditional lithium current sources. 

2. Processes in the materials depending on size of nanoparticles can be accompanied by the 
formation or a two-phase state, or a pseudo-Faraday process with great differential capaci- 
tance. The reason for such a process is quantum-mechanical effect of interference blockade of 
the electronic tunneling into the nonmetallic nanoparticle. However, a mixed mechanism of 
current-generating is not excluded. 

3. Minimization of diffusion resistance in the nanosize silicon dioxide particles promotes the 
achievement of high capacities. 
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4. The phenomenon of current-generating process in nanostructured systems is analyzed on the 
basis of quantum mechanical tunneling. The geometric and energy parameters, that simulate 
the nanostructure, entail the different structure of the electron virtual states. The efficiency 
of tunneling nontrivially depends on the relation between these states and the energy of 
electrons bombarding the nanostructure. Thus, the tunneling efficiency is a non-monotonous 
function on the electron energy. 
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KBaHTOBO-MexaHinHe HBi/ime TyHe/iK)BaHH5i i ecpeKTi/iBHicTb 
cfjapafleeBCbKoro crpyMoyTBopem-m b nopi/iCTHX HaHOCTpyKTypax 

B.A. JlyKiaHei^b, fl.B. MaTy/iKa, I.I. TpnropHaK 

HaL4ioHa/ibHnfi ymBepcuTeT "/IbBiBCbKa ncniTexHiKa", By/i. C. BaHflepw, 12, 79013 JlbBiB, YKpaTHa 

B po6o~ri flocnifl>KeHO TepMOfli/iHaMiMHi Ta KiHeTM h Hi 3aKOHOMipHoc~ri npoL4ecy .nmeBoT iHTepKajismiT C-SiC>2 
HaHOKOM no3i/i~ri b. BcTaHOB/ieHa 3ajie>KHicTb ak p,r\n flwcpepeHL4ii/iHoT eMHOcri, TaK i fljia Ki HeTH h h hx napaMe- 
TpiB npoi4ecy Bifl po3MipiB Ha hokom no3HTiB. 3riflHO 3 pe3y/ibTaTaMH TeopeTHHHoro aHa.ni3y, OTpn Ma Hi 3a.ne- 
>KHOCTi e HacniflKOM KBaHTOBO-Mexa Hi h Horo ecpeKTy iHTepcpepeHLiiCiHoT 6jioKaflw ejieKTpoHHoro TyHejiKDBaHHn 
b HeMeTa^inHy Ha HonacTi/iH Ky. HaBefleHa iMneflaHCHa MOfle/ib npoi;eciB Ta napawieTpi/iHHa ifleHTi/icpiKai4i;R. 
llpeflCTaBJieHa HOBa TexHO/ioria e/ieKTpoxiMinHi/ix reHepaTopiB eHeprfi. 

K/iKJHOBi c/iOBa: h3hoo6'£kt, eneKTpoHHi/itt craH, TyHentoBaHHa, iHTepKannujn 
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